Detectability of large-scale power suppression in the galaxy distribution 
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Suppression in primordial power on the Universe's largest observable scales has been invoked 
as a possible explanation for large-angle observations in the cosmic microwave background, and is 
allowed or predicted by some inflationary models. Here we investigate the extent to which such a 
suppression could be confirmed by the upcoming large- volume redshift surveys. For definiteness, we 
study a simple parametric model of suppression that improves the fit of the vanilla ACDM model 
to the angular correlation function measured by WMAP in cut-sky maps, and at the same time 
improves the fit to the angular power spectrum inferred from the maximum-likelihood analysis 
presented by the WMAP team. We find that the missing power at large scales, favored by WMAP 
observations within the context of this model, will be difficult but not impossible to rule out with 
a galaxy redshift survey with large volume (~ 100 Gpc 3 ). A key requirement for success in ruling 
out power suppression will be having redshifts of most galaxies detected in the imaging survey. 



I. INTRODUCTION 

Measurements of the angular power spectrum of the 
cosmic microwave background (CMB) anisotropics from 
the WMAP experiment have been used to constrain the 
standard cosmological parameters to unprecedented ac- 
curacy [IH1]. At the same time, several anomalies have 
been observed, one of which is the missing power above 
60 degrees on the sky in the maps where the galactic 
plane has been masked [H [5J E] - This is unexpected not 
only because skies with such lack of large-scale power 
are expected with the probability of about 0.03% in the 
standard Gaussian, isotropic model [SHZ], but for two 
other reasons. First, the missing power occurs on the 
largest observable scales, where a cosmological origin is 
arguably most likely. Second, missing correlations are 
inferred from cut-sky (i.e. masked) maps of the CMB, 
which makes the results insensitive to assumptions about 
what lies behind the cut. For review of the missing cor- 
relations (and other so-called "large-angle anomalies" in 
the CMB), see [8]; for debate on this issue, see [81411J: for 
signatures of the anomalies in future polarization obser- 
vations, see [12] . 

In this paper we study the possibility that the primor- 
dial power spectrum is suppressed at large scales. This 
explanation has been invoked before in order to explain 
the low power in the multipole spectrum (e.g. |13j). In 
the meantime, observations have made it apparent that 
the harmonic-space quadrupole and octopole are only 
moderately low (e.g. [21 [15]), and it is really a range 
of low multipoles that conspire to produce the vanishing 
C(9) . Specifically, as discussed in [B] , there is a cancella- 
tion between the combined contributions of C2,---,C§ and 
the contributions of Ce with £ > 6. It is this conspiracy 
that is most disturbing, since it violates the independence 
of the d of different £ that defines statistical isotropy. 

Note however that it is a priori not at all clear that sup- 
pression in the large-scale power can explain the WMAP 
observations on large scales. While the missing large- 
angle correlations in the angular two-point correlation 
function of the CMB C(9) could be trivially explained by 



the missing primordial power, a large suppression would 
lower the harmonic power spectrum Ci, inferred using 
the maximum-likelihood estimator, too much to be con- 
sistent with observations. [We discuss this in Sections IV 
and |v| below.] 

In this paper we perform a two-pronged analysis. First, 
we adopt a simple parametric model for the suppres- 
sion, and perform a detailed analysis to find the sup- 
pressed power spectrum that improves the fit of the 
vanilla ACDM model to the angular correlation function 
measured by WMAP in cut-sky maps, and at the same 
time improves the fit to the angular power spectrum in- 
ferred from the maximum-likelihood analysis presented 
by the WMAP team. Second, we address the following 
question: if the CMB observations are telling us that the 
three-dimensional primordial power spectrum is indeed 
suppressed at large scales (and our adopted model for the 
suppression is at work), could this effect be confirmed in 
redshift surveys, with observations of suppressed cluster- 
ing of galaxies on the largest scales? 

It is important to note that we do not concern our- 
selves with questions recently discussed in the literature 
as to whether the full-sky or the cut-sky measurements 
are more robust. It could be the case that one of these 
measurements, full-sky or cut-sky, is correct while the 
other is not due to some type of systematic error; it could 
also be that both of these measurements are correct (in 
which case the assumption of statistical isotropy is ar- 
guably on less firm footing). We consider these possibili- 
ties separately in order to get a rough idea on what scale 
the data favor power suppression in either case. Regard- 
less of which of these possibilities is true, however, our 
results regarding the detectability of power suppresion, 
presented in Sec. |VI| are valid. 

The paper is organized as follows. In Sec. [TT] we do a 
preliminary investigation in which we attempt to recon- 
struct the suppression of the primordial power spectrum 
(and, correspondingly, the matter power spectrum) di- 
rectly from CMB angular power spectrum measurements 
(the Ci). In Sec. Ill we take the complementary ap- 



proach, parameterizing the suppression and finding how 



it affects the Cg and C{9). Sec. 



IV quantifies how well 
and 



a given suppressed model fits CMB data in both C, 
C(0). Sec. M provides a discussion of the results that we 
obtain from this analysis; we find that large-scale sup- 
pression of power can significantly increase the likelihood 
of the observed CMB at large scales. Finally, in Sec. |VI| 
we discuss the possibility of detecting suppression in the 
matter power spectrum with an upcoming large-volume 
redshift survey. We conclude in Sec. IVII 



II. 



SUPPRESSED POWER: PRELIMINARY 
INVESTIGATIONS 



Let us first review the basic way in which the primor- 
dial power spectrum determines fluctuations in the CMB 
observed today. The CMB temperature anisotropies are 
decomposed into spherical harmonics with coefficients 
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where T is the average temperature of the CMB. The an- 
gular power spectrum, which quantifies the contribution 
to the variance of the temperature fluctuations at each I, 
is then given by the coefficients Cg where, assuming sta- 
tistical isotropy, (ag m ag' m >} = Ci5ai5 mm >. We will also 
consider the angular two-point correlation function 
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where we have assumed statistical isotropy, and the ex- 
pectation is taken over the ensemble of universes. C(9) 
is related to the anisotropy power spectrum by 
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The angular power spectrum Cg is directly related to 
the primordial power spectrum of curvature perturba- 
tions laid down by inflation. The Cg are given in terms 
of the primordial power spectrum by 
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= d(\nk)[Tg(k)] 2 Al(k), 
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where Tg(k) is the transfer function and A%(k) is the 
dimensionless curvature power spectrum 
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where -Pfl(fc) is the curvature power spectrum which, at 
late times and on subhorizon scales, is related to the mat- 
ter density power spectrum P(k) via Pn{k) oc fc~ 4 P(fc). 
We would like to infer the primordial curvature power 
spectrum A 2 R {k) given the angular power spectrum Cg 



measured from the CMB. There are two approaches we 
could take to dealing with Eq. Q: the inverse problem 
(discussed in this section) and the parametric forward 
problem of starting with various power spectra and at- 
tempting to fit the Cg (discussed in the next section and 
pursued in the rest of the paper). 

The first option is to directly calculate A 2 R (k) from 
the measured angular power spectrum Cg. This inverse 
problem, where we know the result of the integration 
but not the integrand, is difficult because the primordial 
power spectrum A R (k) is a three-dimensional quantity 
while the CMB angular power spectrum Cg is a two- 
dimensional, projected quantity. When the problem is 
discretized, as described below, it becomes clear that 
the problem is underdetermined and ill-conditioned, as 
is typical for inverse problems: small changes in the ob- 
served Cg typically lead to large changes in the inferred 
A|(fc). 

Since we are examining the phenomenon of low power 
on large angles in the CMB, it is really the C(8) data 
that we wish to be faithful to, so we take C(9) as our 
starting point rather than Cg. We start from the pixel- 
based measurement of the angular correlation function 
(adopted from [7]), which we denote with a tilde, C(9). 
In order to smooth out the noise in the measured C(9), 
and thereby simplify the inverse problem somewhat, we 
use a "smoothed model" for C{9) that is designed to 
agree with ACDM at small angular scales while closely 
matching the actual WMAP data at larger angular scales. 
To this end, we take the ACDM C(9) and modify it so 
that it smoothly transitions to zero for 9 above roughly 
60 degrees (Figure [l]). 

Inverting Eq. (|3j , we can determine the angular power 
spectrum coefficients Cg inferred from our (smoothed) 
pixel-based estimate C(9) 



Cg = 2ir 



Pg(cos9)C(9)d{cos9). 
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We are now in a position to directly address the inverse 
problem in Eq. (pi). We solve this numerically by dis- 
cretizing the integral 
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The kernel F lk is extracted from CAMB [16]. The basic 
strategy here, distilled in diagrammatic form, is to start 



from C{9) to find the corresponding A|,(fc): 



C{9) -><%-► A 2 R (k). 



(9) 



We attempted two different methods of solving this 
inverse problem, explained more fully in Appendix [A] 
Both methods give similar results, shown with a sam- 
ple reconstruction in Fig. [9] in the Appendix. While 
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FIG. 1. The angular two-point function C(9). Measurement 
from the W-band WMAP 7- year maps, adopted from [7], is 
shown here in blue, while our model in which C{9) smoothly 
transitions to zero at higher 8 is shown in green. This 
smoothed model mimics and idealizes the behavior of the mea- 
sured C{6) at large scales but follows the ACDM prediction 
(red curve) on smaller scales, below roughly 50 degrees. The 
dotted line shows zero correlation for reference. 



this result can only be suggestive (given that we used 
a smoothed model for C(9) and given that the inverse 
problem is ill-conditioned and unclcrdctcrmined), it does 
indicate that a transition to low/zero power on large an- 
gular scales in C{&) can be explained by suppression at 
low k in the primordial power spectrum A^,(fc). If the 
transition to zero power in C(8) occurs at about 60 de- 
grees, as it appears to do in the WMAP cut-sky data, 
then this corresponds to power suppression at scales of 
k < 10- 3 - 6 Mpc _1 « 3.5 x lCT 4 /i/Mpc. 



III. SUPPRESSED PRIMORDIAL 
LARGE-SCALE POWER 

The inverse approach from the previous section and 
the Appendix shows that the direct inversion of our 
smoothed C(9) leads to a suppression of P(k) at 
log 10 (fc c /(/i/Mpc)) < —3.5, but the inversion is very 
noisy and non-robust, as expected. We now change tac- 
tics and move to the alternative approach to Eq. Q. 
Instead of treating this as an inverse problem, we now pa- 
rameterize An(fc) and treat this as a (much more stable) 
forward problem. We utilize a three-parameter model pa- 
rameterizing the suppressed A R (k) with an exponential 
cutoff: following [13] and [17], we write 
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(10) 
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where we have implicitly defined the factor 
S{k;k c ,a,(3) = 1 - /3cxp(-(fc/fc c ) Q ) by which the 
power spectrum is suppressed. The parameter k c 
controls the fc-value of the transition; a controls the 
sharpness of the transition; and the extra parameter 



1. 

S°- 

^ 0. 



log(* c )=-3.35 
log(* c )=-2.85 



0.0 



lo Sio(- 



k 



h/Mpc 



^ o 



t = 1.0 

« = 3.0 



lo §io(- 



h/Mpc 



1.0 






,3 = 0.5 

■ ■■■i ,3 = 1.0 


*v 


?5 








"*01 


: 




^0.2 
0.0 


c 




,.-* 




6 -5 -4 -3 


2 


-1 



lo Sio(- 



/,• 



h/Mpc 



FIG. 2. Illustration of what the suppression factor S(k) looks 
like for various combinations of parameters. The default pa- 
rameters are log(fc c ) = log 10 (fc c /(/i/Mpc)) = —2.85, a = 3.0, 
and P = 1.0. In each plot, one of these parameters is varied, 
while the other two are held fixed at their default values. 



/3, which is not found in [13] or [17], allows the power 
spectrum to plateau to a value other than zero at low 
k (S(k) — > 1 — j3 for k — > 0). Note that this parameter- 
ization has enough freedom to mimic the results of the 
inversion shown in Fig. [9] almost perfectly. 

For a given set of parameters {fc c ,a,/3}, we have a 
well-defined A|j(fc) and can thus use the numerical kernel 
Ffk to find the corresponding C/s, as in Eq. ^. We 
can then determine which combinations of parameters 
give CYs that fit the observed WMAP C/s, and likewise 
the observed C(9). We are now moving in a direction 
opposite the one in Eq. ^: 



A 2 R (k) -^C e ^ C{6). 



(12) 



Plots of the suppression factor S(k) for several sample 
parameter values are shown in Figure [2] 

When varying the suppression parameter k c (and, op- 
tionally, a and /3), we have not simultaneously varied 



other parameters that describe the primordial power 
spectrum, such as the dark matter and baryon densi- 
ties, spectral index, etc. The reason, in addition to sim- 
plicity, is that none of these other parameters can mimic 
the large-scale suppression of power, and therefore, power 
suppression is not degenerate with other cosmological pa- 
rameters. The one possible exception would be primor- 
dial nongaussianity of the local type, which does indeed 
affect the power spectrum of halos (and, thus, galaxies) 
on large scales [18]; however, including this degeneracy 
is beyond the scope of this project. 



IV. STATISTICAL TESTS 

We are interested in how a suppressed primordial 
power spectrum, as given in Eq. (11), affects both Ct 
and C{9). Figures 3] and [I] show examples of how these 
quantities vary with changes in the parameters. 

More specifically, we look to quantify whether, and 
to what extent, a suppressed primordial power spec- 
trum gives a better fit to observations of Cg (typically 
inferred using maximum-likelihood-type techniques at 
large scales) and C(9) estimated on cut-sky maps using 
a pixel-based estimator. We restrict attention here to 
varying k c rather than all three parameters, since vari- 
ations in k c have the greatest effects on the likelihood, 
and also since k c directly controls the scale at which the 
suppression occurs. 

In the case of Cg, it is relatively straightforward to 
quantify the fit between a given suppressed model and the 
WMAP data. We simply generate a suppressed model 
using CAMB and then feed the resulting Cg spectrum 
into the WMAP likelihood codqj to obtain a goodncss- 
of-fit criterion. Further details are covered in Sec. IIV Al 

Quantifying the fit between a suppressed-model C(0) 
and the WMAP data requires a bit more work. In Sec. 



IV B we examine the statistic Si/%, which gives a mea- 
sure of correlations in the CMB above scales of 60°. In 
analogy to our treatment of the angular power spectrum, 
we define a \ 2 statistic to quantify the goodness of fit be- 
tween a given suppressed model and the WMAP data. 
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FIG. 3. Effects of the suppression on the angular power 
spectrum Cg. The default parameters are again log(fc c ) = 
log 10 {k c /(h/Mpc)) = -2.85, a = 3.0, and f) = 1.0. In each 
panel, one of these parameters is varied, while the other two 
are held fixed at their default values. The WMAP measure- 
ments are shown as points with measurement error bars (too 
small to see in this plot). Cosmic variance (assuming full-sky 
measurements) is plotted as a band around the most heavily 
suppressed model (green curve) in each panel. 



A. Angular Power Spectrum Ct 

To quantify how well the suppressed angular power 
spectrum Cg ySup fits the WMAP observations Cc,wmap, 
we use the WMAP likelihood code to compute 
Xc e ,su P ( k c) = -21n£ sup (fc c ) from Cg ySup . We likewise 
compute Xc e ,unsuv = — 21n£ unsup from the unsuppressed 
ACDM power spectrum CV jlmsup , finding the difference 



1 Available from http://lambda.gsfc.nasa.gov/product/map/ 
dr4/likelihood_info.cfm 



A Xc t ( k c) = X 2 C t , S up( k c) - Xc f _ 



,unsup 



(13) 



as a final quantification of how well the suppressed model 
fits Cg data relative to the unsuppressed ACDM model. 

A negative A%^ indicates that the suppressed model 
is a better fit to the WMAP data than ACDM . For mod- 
els with a large amount of suppression (i.e., high k c ), the 
suppressed model affects the Cg at increasingly smaller- 
scales (higher £), making them inconsistent with the Cg 
measured from WMAP and thus making the quantity 
A^c" large and positive. 

The results are shown further below in Fig. [6l with 
discussion in Section fVl 
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FIG. 4. Effects of the suppression on the real-space two- 
point correlation function in the CMB, C(8). The default 
parameters are again log(fc c ) = log 10 (fe c /(/i/Mpc)) = —2.85, 
a — 3.0, and /3 — 1.0. In each plot, one of these parameters 
is varied, while the other two are held fixed at their default 
values. 



B. Angular Correlation Function C(d) 

We now study to what extent the power suppression 
from our model can account for the low correlations ob- 
served above 60° on the cut-sky WMAP maps [U [6] . 
The statistic Si/ 2 , defined in pQ, quantifies the lack of 
correlation above 60 degrees 



Si 
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It is possible to calculate S\/2 directly from the Cn 's: 
S1/2 = ^^2(2e+l)(2£' + l)Ceh,i'(l/2)C e/ . (15) 



W 



For details of how the quantity I^gj is calculated, see 
Appendix A of the published version of Ref . [5] . 
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FIG. 5. The distribution of our realizations of the statistic 
Si/2 for the unsuppressed (red histograms) and a sample sup- 
pressed (green histograms) model. The suppressed model has 
logio(fc c /(V M Pc)) = -2.85, a = 3.0, and = 1.0, and has 
a better fit than the unsuppressed model by Axl, ,„ = —7.9. 

1/2 

The bottom panel clearly shows that the distribution of Si/ 2 
is lognormal, whether or not the underlying power spectrum 
is suppressed. 



We would like to generate statistical realizations of the 
angular power spectrum based on the underlying primor- 
dial power spectrum A|,. To do that, we first calculate 
the expected angular power spectrum C^ )0r i g inal (using 
Eq. (nl), and then create realizations using 
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where the multiplicative factor 

r(fc = (2i + i)/ sky /2, 

(2^+l)/ sky 
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= 2) 



(16) 



(17) 



the numerator is drawn from a gamma distribution with 
scale parameter 9 = 2 and shape parameter k = (21 + 
l)/sky/2, and the denominator ensures that the mean of 
/ is unity. The reason we draw from a gamma distribu- 
tion is that this is the appropriate generalization of a \ 2 
distribution when the number of degrees of freedom is 
non-integer, as it is above for / s k y 7^ 1; T(k = r/2, = 2) 
is identical to a \ 2 distribution for integer r degrees of 
freedom. We adopt / s k y = 0.75 for the remainder of this 
paper. [Modeling of the noise is unnecessary since cosmic 
variance dominates at these large scales.] 

We examine the resulting distribution of S1/2 values by 
performing 100,000 realizations of the Ct for 2 < t < 50 
(going to higher values of £ barely changes the Si/ 2 statis- 
tic, since scales above 60° are mostly affected by i < 10), 



assuming central values C£,originai calculated based on the 
suppressed primordial power spectrum as in Eq. (I8J) , and 
calculating Si/ 2 for each set of (^realization- We find that 
S1/2 is distributed approximately according to a lognor- 
mal distribution, both for suppressed and unsuppressed 
models, and regardless of the particular value of k c . This 
is illustrated in Fig. [5] 

In the sample suppressed model (log 10 (fc c /(/i/Mpc)) = 
—2.85) shown in Fig. [5| the histogram of Si/ 2 peaks at 
1000([iK ) 4 , has a mean of 8300(/iAT) 4 , and a median of 
4300(/iAT) 4 . The histogram of m(S 1/2 ) peaks at 8.4, cor- 
responding to an S1/2 of 4400(^AT) 4 . These values are all 
much lower than the mean value expected in the best- 
fit ACDM cosmology (about 50,000(iiK) 4 ), but bigger 
than the value measured in WMAP cut-sky maps (about 
1000(^?f) 4 ). 

In order to calculate a x 2 statistic in analogy to the 
Axq above, we first transform the lognormal Sx/ 2 dis- 
tribution to a Gaussian by taking the natural log of the 
£1/2 values. The result, a nearly perfect Gaussian, is 
shown for the unsuppressed and the sample suppressed 
model in the lower panel of Fig. [5j We can then calcu- 
late the x 2 corresponding to the probability of getting a 



certain value S°^ of S\/ 2 
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where ln(Si/2{k c )) is the mean over the realizations of 
hi(Si/ 2 ) for the given k c and <Ti n (s 1/2 ) is the standard de- 
viation over all realizations. For the purposes of this pa- 
per, we choose S°^ = 1000 (^K) 4 , since this is (roughly) 
the value of Si/ 2 favored by the cut-sky WMAP obser- 
vations [IJI5H7]. 

We also performed 6,500,000 realizations of the Ce as- 
suming central values C^originai corresponding to the un- 
suppressed ACDM model. From the S1/2 values associ- 
ated with these realizations, we calculate y% ,„,„„„ in 

"-'l/2i UIlSlip 



exact analogy to Eq. (18), and then compute 



A xi 1/2 (fcc) = x! 1/2 ,su P Oc) - Xs 1 



1/2 > unsu P' 



(19) 



A combined statistic that takes into account both the 
measurements of the angular power spectrum Ce and the 
total angular correlation above 60° parameterized by the 
statistic 5*1/2, is then given bjr] 



2 Since we are interested in how likely the low value of Sj 12 ~ 
1000((j,K) 4 is, given suppression of power, we could have simply 
calculated P(Si 12 < 1000) - the probability that Si 12 is as low 
as 1000 - instead of performing the more complicated calculation 
above to obtain i- > (5 1 / 2 )- However, the danger in doing this 
is that suppression on small enough scales leads to values of 
S1/2 that are much lower than 1000, and then the probability 
that S1/2 is as high as 1000 should become low. Considering 
the Gaussian likelihood in ln(S]/2), as we have done, correctly 
penalizes values of Si 12 that are too low or too high. 
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Both P(Si/ 2 ) and P(Ce) are normalized so that their 
values for the unsuppressed ACDM model are 1. Hence 
P(Si/ 2 ) and P{Cg) should be interpreted as the improve- 
ment (relative to fiducial unsuppressed ACDM) in how 
well a given suppressed model fits the WMAP data for 
Ce and S\/ 2 . Note that we are not taking the correlation 
between the (maximum-likelihood) Ce and the (pixel- 
based) S1/2 into account in our statistic C{k c ). We define 
the statistic in the simplest possible way, by multiplying 
the individual likelihoods in Ce and S\/ 2 . This simple 
combination is sufficient, since it favors suppression on 
scales between the scales which P(Si/ 2 ) and P(Ce) inde- 
pendently prefer (this is confirmed in Fig. pi), and thus 
captures the essence of how these two quantities jointly 
favor suppression. Note that the main results of this pa- 



per, presented in Sec. VI do not depend on exactly how 



we combine the likelihood in the measured full-sky Ce 
and cut-sky Si/ 2 ; we only use C{k c ) from Eq. p0|) to 



get a rough idea about the scale at which suppression is 
favored by the data if both measurements are taken at 
face value. We then proceed in Sec. |Vl]to study the de- 
tectability of suppressed power corresponding to a range 
of values of the suppression scale k c ; these results do not 
depend on how P{Ce) and P{Si/ 2 ) are combined. 

We are now in a position to determine which values of 
the suppression scale k c improve the joint fit to the angu- 
lar power spectrum Ce and the cut-sky measurements of 
the angular correlation function C{9) above 60° (quanti- 
fied by the statistic Si/ 2 ). 



V. CURRENT CONSTRAINTS FROM THE 
CMB 

The top panel of Fig. [6| shows P{Si/ 2 ), P(Ce), and 
their product as a function of k c , with a held fixed at 
3.0 and j3 held fixed at 1.0. The bottom panel displays 
the same result using Xs an d Xct on the vertical axis 
instead of P(S 1/2 ) and P{C e ). 

As indicated in Fig. [6j introducing suppression in 
the primordial power spectrum can increase the likeli- 
hood of both the observed Cg and the observed Si/ 2 , 
but these two observations favor suppression at differ- 
ent scales. The likelihood of the Ce is improved by at 
most a factor of 2.2, with the improvement peaking at 
log 10 (fc c /(/i/Mpc)) = —3.4, while greater suppressions 
can improve the likelihood of the S , 1 / 2 data by huge fac- 
tors of up to 131, peaking around log 10 (A; c /(/i/Mpc)) = 
—2.6 (note the plotting scale of likelihoods in Fig. [6j 
where individual likelihood curves are divided or mul- 
tiplied by 10 for visual clarity). The Ce measurements 
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FIG. 6. The improvements in P(Si/2), P(Ce), and the prod- 
uct thereof - all relative to the unsuppressed model - are 
plotted as a function of k c . The top panel shows P{Si/z) 
divided by 10 and P(C t ) multiplied by 10. The P(S 1/2 ) 
values are based on 100,000 realizations of the Ce's as de- 
scribed in the text. The bottom panel shows the same infor- 
mation as the top panel, but puts it in terms of chi-square 
goodness-of-fit statistics. The maximum improvement in 
P (S'1/2) x P(Ce), relative to the unsuppressed ACDM model, 
is a factor of 8.7 (Ax| + ^Xc, — — 4.3), occurring at 
log 10 (fcc/(/i/Mpc)) ~ —3.3. All calculations were performed 
assuming / s k y = 0.75. 



thus favor suppression on very large scales, while the cut- 
sky S1/2 favor suppression all the way down to relatively 
small scales, where suppression is overwhelmingly ruled 
out by Ce data. This is another reminder of just how low 
the pixel-based cut-sky measurement of S1/2 is. It is also 
a reminder of the fact that such a low value of S'1/2 rep- 
resents a conspiracy of the low-£ Ct values: the WMAP 
cut-sky data indicate that S'1/2 is sufficiently low as to 
strongly (by factors of over 100 in likelihood) favor sup- 
pression of primordial power at scales corresponding to 
k = 10~ 2 ' 6 « 0.003 /i/Mpc, even though the maximum- 
likelihood Ct favor suppression only weakly, at far larger 
scales, and overwhelmingly reject the possibility of sup- 
pression at the scales favored by the cut-sky S'1/2. A sky 
with such a low Si/ 2 as the WMAP cut sky ought to have 
CiS that are even more suppressed than the most sup- 
pressed model in Fig. [3] What we see instead are \aw-(. 
Ct values that are not so close to zero, but which instead 
conspire with one another in just such a manner as to 
produce an exceptionally low value of Si/ 2 anyway [5J. 
In any case, we have calculated the product statistic 



as a measure of how well a given suppressed model fits 
the WMAP data in both Ct and C(8) (the latter via 
the specific statistic Sj/2). Since Ct and S1/2 data favor 
suppression at such different scales, there should be a 
"sweet spot" somewhere between the peak in P{Ct) and 
the peak in P(Si/ 2 ), where suppression is moderately fa- 
vored by both Ct and S'1/2, or heavily favored by one and 
still allowed by the other. This is indeed what we find, as 
indicated by the red curve in Fig. [6] Because suppression 
on overly small scales (below log 10 (A; c /(/i/Mpc)) ~ —3.2) 
brings the Ct data for the suppressed model into se- 
vere conflict with the WMAP CVs, the peak of the £ 
curve occurs above these scales, even in spite of the huge 
gains in likelihood that P{S\/2) gives us at much smaller- 
scales, where the gain in P(Si/ 2 ) is still substantial and 
the Ct data still favor - or at least do not heavily dis- 
favor - suppression. The maximum improvement pos- 
sible in P(Si/ 2 ) x P(Ct), relative to the unsuppressed 
ACDM model, is a factor of 8.7 (A x 2 Sl 2 +&X 2 Ce = ~ 4 - 3 )' 
occurring at log 10 (A; c /(/i/Mpc)) ~ —3.3. 

The WMAP likelihood code uses a Bayesian (Gibbs 
sampler) maximum likelihood method (e.g. |19j ) to com- 
pute the fiducial Ct's at the multipoles £ < 32 (5U1 |2"T] . 
We experimented with running the likelihood code using 
pseudo-CV estimates at low multipoleaj and discovered 
that in this case, suppression is much more heavily fa- 
vored by the Ct likelihood than it is in the (presumably 
more accurate) Gibbs sampler, or else a similar Maxi- 
mum Likelihood Estimate (MLE) method. This result 
is expected, and holds because the Ct's that result from 
the pseudo-C^ estimates are lower than those found us- 
ing the Gibbs sampler method (see e.g. Fig. 15 in [15)). 
and suppression fits them better. In this case we can get 
AXc as l° w as — 7.6, corresponding to improvements in 
P(Ct) by factors of up to 44 (as opposed to roughly 2 in 
the best-case scenarios discussed above). 



VI. FUTURE DETECTABILITY USING 
GALAXY SURVEYS 



The results of the previous section indicate that sup- 
pression of primordial power on large scales can increase 
the likelihood of both the observed Ct angular power 
spectrum and the observed cut-sky value of S'1/2, pro- 
vided the suppression "kicks in" on appropriate scales. 
Now we turn to the question of whether large-scale sup- 
pressed power could be detected in the matter power 
spectrum as measured by upcoming redshift surveys such 
as the Large Synoptic Survey Telescope (LSST; [H]). If 
the zero-correlation signature of large-angle C(9) in the 
CMB is an authentic effect indicating a deficit of power 



3 We did this by turning off the use_lowl_TT option in the test .F90 
routine of the WMAP likelihood code, and also switching off 



C(k c ) = P{S 1 / 2 ) X P{Ct) , Or alternatively A%| +AXc, > polarization by turning off the use.TE and use_lowl_pol options. 



on the Universe's largest scales, is it possible to cross- 
check and verify this result using large-scale-structure 
data? 

Given suppression of the prim ord ial power spectrum 
A|j(fc) as parameterized in Eq. (Ill, the matter power 
spectrum will be suppressed by the same factor as A|.(fc): 



P 3up(Kj — i>(fc)-r unsU p(fC) 



(21) 



where S(k) = S(k; k c , a, f3) is the same as before. We 
wish to determine whether this suppressed matter power 
spectrum could be distinguished from the unsuppressed 
ACDM matter power spectrum P un sup(&) by a large- 
volume redshift survey. 

When measuring the matter power spectrum with a 
redshift survey, the error bars in each thin slice in redshift 
dz and wavenumber dk are given by the Feldman-Kaiser- 
Peacock (FKP; [35]) formula 



ap(k,z) 



4ir 2 P(k,z) 2 
k 2 dkdV cS 



(22) 



where the effective volume element dV c s(k,z) is related 
to a comoving volume element via 



dV eS (k,z) = 



n(z)P(k,z) 
l + n(z)P(k,z) 



dV s , 



,{z). (23) 



The differential survey volume is given in terms of dz via 



r(z) 2 

^'survey — ^'survey «-/ ~(*>^i 



H(z) 



(24) 



where r(z) is the comoving distance as a function of red- 
shift, H(z) is the Hubble parameter as a function of red- 
shift, and f2 S urvcy is the angular size of the survey in 
steradians. 

The number density of galaxies n(z) can be found from 



n(z) = m{z) x 



N tl 



n 



survey 



fm(z)[r(z) 2 /H(z)}dz 



(25) 



where the second term on the right-hand side provides a 
normalization. Here iV-tot is the total number of galaxies 
in the survey and m(z) is the (unnormalized) number 
density of galaxies, whose functional form we adopt to 
be 



Given a suppressed power spectrum as in Eq. (21 ), we 
can calculate 



d\ 



2 [-* unsupv™'; %) ■*sup\"') %)\ 



(27) 



nP 

1 + nP 



O k 3 

2 ^survey "* 



4tt 2 P 2 



x\P 



unsup 



-Rmp] d(\nk)dz 



(28) 



and then integrate in order to find the x 2 statistic for how 
well the survey can distinguish between the suppressed 
and unsuppressed models. 

Note that in the above two equations, wherever a 
P = P(k, z) occurs without being marked as either sup- 
pressed or unsuppressed, this is intended to indicate that 
either P sup or P unS up m ay be used. Whether we use P sup 
or Punsup depends entirely on which question we are try- 
ing to answer: If we use suppressed-model error bars, 
then this % 2 = x 2 up indicates at what confidence level 
the survey can rule out suppression. Meanwhile, if we 
use unsuppressed- model error bars, then x 2 = xSnsup m ~ 
dicates at what confidence level the survey can rule out 
the unsuppressed ACDM model. Ruling out ACDM is 
considerably more ambitious than ruling out suppres- 
sion, since the error bars tend to be smaller when they 
are based on the suppressed model (due to the fact that 
op oc P + 1/n). 

This is illustrated in Fig. [7J The plot shows the un- 
suppressed matter power spectrum, along with the sup- 
pressed version for a particular choice of parameters. The 
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(26) 



We take zq = 0.35, corresponding to the density roughly 
expected in the imaging portion of the LSST survey 
[24], and assume a 23,000 square degree redshift survey 
with 0.5, 5 or 50 spectra per square arcminute. [Note 
that the 0.5 and 5 gal/arcmin 2 cases are realistic, be- 
ing targeted by surveys in the near future [25] [26] , while 
50 gal/arcmin 2 corresponds to the more aggressive case 
where spectra of most galaxies in the imaging portion of 
the survey are taken.] 



FIG. 7. Matter power spectrum P(k) with and without 
suppression. In the suppressed model, the parameters are 
logi (fc c /(V M Pc)) = -2.85, a = 3.0, and f3 = 1.0. The 
power spectrum is shown at z = for a redshift survey with 
50 galaxies per square arcminute. The error bars are based 
on the suppressed power spectrum, which means that a suf- 
ficiently high \ 2 value here would indicate the possibility of 
ruling out suppression. The value of \ 2 turns out to be 57.8, 
good enough to rule out suppression at roughly 7.6a. (Bins 
without vertical error bars contribute nothing to the \ 2 ~ i n 
effect, their error bars are infinite.) Note the log scale on both 
axes. 
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FIG. 8. The detectability of suppression as a function of k c . 
These results apply to a survey that extends from z = to 
z = 3 covering 23,000 square degrees of sky. From bottom to 
top in each set of lines, we assume 0.5, 5 or 50 (spectroscopic) 
galaxies per square arcminute. [Note that the 0.5 gal/arcmin 2 
case is entirely realistic in the near future, corresponding to 
the number density of spectra planned by e.g. BigBoss [25] . 
while 50 gal/arcmin 2 corresponds to the more aggressive case 
where spectra of most galaxies in a large- volume imaging sur- 
vey are taken.] Here a is fixed at 3.0 and /3 is fixed at 1.0. 



goal of calculating \ as m Eq. (28) is to determine 



whether the unsuppressed power spectrum can be distin- 
guished from the suppressed power spectrum for a given 
set of parameters within the error bars that would be 
set by an LSST-like survey. For the case pictured - in 
which the unsuppressed model is taken as true, and the 
error bars are calculated based on the suppressed model 
which is being tested - it is possible to rule out suppres- 
sion with high statistical significance. The opposite is, 
however, not true: if the suppressed model is true, it will 
be very difficult to rule out the standard unsuppressed 
ACDM due to its larger errors. 

For example, survey measurements that fell along the 
curve predicted for the unsuppressed P(k) would, in the 
case shown in Fig. [7J fall outside some of the suppressed- 
model error bars, and ultimately combine to give a to- 
tal x 2 °f 57.8, given the survey parameters outlined in 
the next paragraph. Meanwhile, taking the unsuppressed 
model as fiducial would allow for the possibility of survey 
measurements ruling out ACDM, but the total x 2 would 
shrink to 2.0 due to the larger error bars. 

The final results for the detectability of suppression 
are shown in Fig. [8] Instead of plotting \ 2 we show the 
number of sigmas (i.e. y'Ax 2 ) at which suppressed and 
unsuppressed power spectra can be distinguished assum- 
ing one degree of freedom on the measurements of P(k). 
The figure shows the results as a function of k c , holding 
the parameters a and /3 fixed at 3.0 and 1.0, respectively, 
and assuming three different possible values for the num- 
ber of galaxies observed per square arcminute (0.5, 5, 
and 50) in the spectroscopic survey. We also examined 
the results with different values of a and /3, but changes 
in these parameters do not greatly affect the results un- 



less /3 becomes close to zero. The scale of the suppression 
as determined by k c is by far the greatest contributing 
factor in determining whether a given suppressed model 
will be detectable to a large-volume redshift survey. 

Comparison of Fig. [8] with Fig. [6] shows that if present 
data for Ci and C(9) truly point to suppression of the 
primordial power spectrum, that suppression is likely on 
scales that are too large for foreseeable redshift surveys 
to either detect or rule out. The most optimistic scenario 
shown in Fig. [HJ in which there are 50 galaxies per square 
arcminute in the spectroscopic survey, still cannot (at 
3er) rule out suppression if log 10 (fc c /(/i/Mpc)) < —3.0, 
and cannot rule out ACDM unless the Universe actu- 
ally shows suppression of the matter power spectrum on 
much smaller scales, with log 10 (fc c /(/i/Mpc)) > —2.7. 
Suppression on scales this small is strongly disfavored 
by WMAP Ce observations. Meanwhile, the scales on 
which WMAP observations tend to favor suppression 
(log 10 (fc c /(/i/Mpc)) ~ —3.3) are nearly inaccessible to 
galaxy surveys. This is a reflection of the fact that the 
CMB probes much larger scales than even the largest- 
volume redshift surveys of the near future. 

If only the cut-sky Si /2 statistic is taken into account, 
CMB observations heavily favor suppression on scales 
where suppression would be readily detectable by red- 
shift surveys, at several sigma, for number densities of 
galaxies expected in near-future spectroscopic samples. 



VII. CONCLUSIONS 

In this paper we have studied the suppression of pri- 
mordial power on large scales as a possible explanation 
for the CMB observations. Without considering partic- 
ular physical models for the suppression, we adopted a 
more pragnratic approach and addressed the following 
question: do the suppressed models actually improve the 
likelihood of the observed CMB sky and, if so, can the 
upcoming large-volume galaxy redshift surveys be used 
to confirm this suppression? 

We first motivated our search by attempting to invert 
the observations of the angular power spectrum Cg, in or- 
der to reconstruct the three-dimensional power spectrum 
P(k). As expected, this procedure is very unreliable and 
noisy due to the nature of the inverse problem; never- 
theless, we obtained useful hints for the form of the sup- 
pression that we should be considering (see Fig. [9] in the 
Appendix) . 

We then proceeded to use a parametric model of the 



suppression (Eq. (11)), with the most important param- 
eter (and the only one we varied in our analysis) being 
the suppression scale k c . We found (see Fig. [6| that the 
angular power spectrum Cg, traditionally inferred using 
maximum-likclihood-type estimators, prefers a moderate 
suppression of power; conversely, the cut-sky pixel-based 
correlation C(0) prefers a stronger suppression. It is also 
possible that both the full-sky measurement of Cg and 
the cut-sky measurement of C{9) are not anomalous, but 
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rather that the underlying cosmological model is not sta- 
tistically isotropic. While it is not clear how to write 
down the combined likelihood in the full-sky and cut- 
sky measurements without assuming statistical isotropy, 
our simple choice (Eq. (20)) prefers the suppression at 
log 10 (fc c /(/i/Mpc)) ~ —3.3, and increases the combined 
likelihood by about a factor of 8.7, corresponding to 
Ax 2 = -4.3. 

Detectability of such a large-scale suppression with 
future surveys will be difficult, however, as shown in 
Fig. [8] In order to detect the suppression favored by 
the CMB angular power spectrum, an LSST-type sur- 
vey, with a volume of about 100 Gpc 3 and a very large 
number of galaxy redshifts measured, will be necessary. 
Roughly speaking, a statistically significant ruling-out of 
the power suppression will require spectra taken of most 
galaxies in the imaging portion of the survey; this will 
require a ~ 10-meter ground-based, or a ~ 1.5-meter 
space-based, telescope dedicated to taking spectra. Al- 
ternatively, photometric redshift techniques may some- 
day become so accurate that our preferred case of "nearly 
all galaxies being spectroscopic" is validated relatively 
straightforwardly. 

Additionally, we point out that suppressed power will 
be more easily ruled out (given that the true power is 
not suppressed) than vice versa, essentially because the 
model being tested has smaller cosmic-variance errors if 
it has lower power. Therefore, if indeed we live in the 
universe with the true power spectrum of density fluc- 
tuations being standard inflationary power-law (i.e. un- 
suppressed), then, for example, a survey covering half 
the sky with 5 galaxy redshifts per square arcminute will 
be able to rule out power suppressions on scales above 
roughly 1 Gpc at 3ct confidence; suppression extending 
to smaller scales is even easier to detect. 

Overall, we are optimistic about the prospects of 
galaxy surveys to test models of the suppressed large- 
scale power of primordial fluctuations. Dark Energy Sur- 
vey (DES; [57)), Baryon Oscillation Spectroscopic Sur- 
vey (BOSS; [28]) and, especially, very-large- volume sur- 
veys such as the LSST [22], Joint Dark Energy Mission 
(JDEM; [29]), Euclid [30], and BigBoss [25], will be able 
to test, at least in part, observations of CMB experiments 
on the largest observable scales. 



Appendix A: Direct inversion to obtain the 
primordial power spectrum 

As pointed out in Sec. [TT], the CMB angular power 
spectrum Ci is given in terms of the primordial power 
spectrum by 

l{i+ 2 l )Cl = J d (^ k) [T,(k)] 2 Al(k) (Al) 
= Y,F ik A 2 R (k), (A2) 



where the discretized numerical kernel F^ is extracted 
from CAMB [16 . Trying to find A 2 R (k) from a given set 
of CVs (which themselves correspond to a given C(0)) 
is an inverse problem, which we attempted to solve us- 
ing two different strategies. [We also attempted a third, 
doing a simple matrix inversion of the kernel, but this 
strategy simply does not work due to the extreme ill- 
conditioning.] The first is the Richardson-Lucy method, 
an algorithm that iteratively solves for the portion of the 
sum that multiplies the kernel [3"lT(3"4] : 



A 2 +1 (fc) = A 2 (fc) £ F lk C * 



obs 



ci ' 



(A3) 



where Cg bs are the observed CVs (in this case, the C/s 
correspondingto the smoothed C(8) shown as t he green 
curve in Fig. [I ) , C\ is calculated from Eq. ( A2 1 for each 
iteration i, anc 



F(k — Ftk/ / „ Ft, 



(A4) 
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The method converges to a solution for A„(fc), but has no 
special properties guaranteeing convergence or smooth- 
ness of the solution. 

A second strategy makes use of linear regularization, 
which is one way of putting extra constraints on the 
solution. The angular power spectrum Cg, is a two- 
dimensional quantity, while the primordial power spec- 
trum A 2 (fc) is three-dimensional, and so finding the latter 
from the former is an underdetermined problem. Linear 
regularization compensates for the fact that the inverse 
problem is underdetermined. 

Numerical Recipes [35] outlines one method of regu- 
larizing, which favors a constant solution and penalizes 
deviations from this. The goal here is to solve the equa- 
tion 
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Ku 



(A5) 



where we refer to matrix K as the kernel; b is a known 
vector; and u is the vector to be solved for. Regulariza- 
tion does this by solving the regularized equation 



(K T K + \H)u = K T b, 



(A6) 
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where H is a matrix that takes a different form depending 
on whether the regularization is linear (penalizing devi- 
ations from a constant solution), quadratic, etc., and A 
is a parameter that controls how strong the regularizing 
constraint is: higher values of A impose stronger regular- 
izing constraints on the solution. In our case, we take the 
kernel K to be i 7 ^, u to be a column vector correspond- 
ing to A^j(fc), and b to be a column vector corresponding 
to C° bs . We adopt H corresponding to linear regulariza- 
tion. 
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FIG. 9. Sample results of the direct inversion. The green 
points represent the factor S(k) (see definition of S(k) in 
Eq. (Ill) by which the primordial power spectrum Afj is 
suppressed as determined by the regularized inversion from 
Eq. (A2|. The angular power spectrum Ce used in this inver- 



sion corresponds to the C(9) shown as the smoothed green 
curve in Fig. [I] though we could have in principle used the 
actual WM AP data shown as a blue curve in the same figure. 



It is difficult to get consistent results from either of 
these strategies, given the ill-conditioned and underde- 
termined nature of the inverse problem. This is part of 
the reason why we chose to focus most of our attention 
on doing the forward problem outlined in Eq. ( 12 ). How- 



ever, to the extent that consistent results are possible, 
both strategies give similar solutions. A sample result 
for the suppression factor S(k) (see definition of S(k) in 
Eq. ( 11 )) is shown in Fig. [9J 

The most notable feature of the inversion result is 
that it transitions to near-zero power at large scales/low 
k, with a form suggesting an exponential cutoff; this 
provides motivation for adopting the form we d id for 
parametrizing A|j in the forward problem (Eq.lllj)). Un- 
fortunately, a direct inversion of the sort described in this 
Appendix requires some fine-tuning in order to get results 



of this quality, which is why we have emphasized that 
these results are suggestive rather than conclusive. The 
results in Fig. M (which correspond directly to the regu- 
larized inversion, but are also similar to the results of the 
Richardson-Lucy method) rely on careful tuning of the 
regularization parameter A to ensure that the result does 
not become negative at low k. [If the regularization is 
"not strong enough," with A too small, deviations from a 
constant solution are not sufficiently penalized to prevent 
the solution becoming negative. If the regularization is 
"too strong," with A too large, the solution simply stays 
constant at roughly 1. Only in particular intermediate 
cases does it transition nicely from to 1.] 

In addition to this issue, the fact that the inversion 
results show S(k) of exactly 1 at high A; is a result of 
a mechanism that was put into the solution process by 
hand. Without enforcing the high-fc value, the results of- 
ten converge to a constant which deviates from unity at 
high k. To compensate for this and for the issue that the 
solution does not always asymptote to nonnegative values 
at low k, we attempted a modification of the regularized 
inversion in which deviations from and 1 are penal- 
ized (rather than deviations from a constant solution, as 
in linear regularization). The solutions we obtained us- 
ing this method have the same general form as shown in 
Fig. [9l with a transition from to 1 somewhere between 
log 10 (fc c /(/i/Mpc)) = —3 and —4, but the results are even 
noisier than the results obtained with Richardson-Lucy 
and linear regularization. 

When the WMAP data was used directly as input, 
rather than the Ct corresponding to the smoothed model 
in Fig. [l] solutions to the inverse problem were even nois- 
ier. 

Finally, neither the regularized inversion nor the 
Richardson-Lucy method give error bars with which the 
precision of the inversion might be judged. For all these 
reasons, the results of the direct inversion cannot be 
taken as anything more than suggestive. With that 
caveat, it is still notable that results of both the regu- 
larized inversion and Richardson-Lucy method do con- 
sistently suggest a transition from suppressed power at 
low k to unsuppressed power at high k. This provides a 
hint of the fact that the likelihood of the WMAP Cg, and 
C(6) data may be increased by introducing suppression, 
as explored much more fully, and confirmed, in Sections 
En and M 
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